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Abstract 

We introduce a novel framework for the stability analysis of discrete-time linear switching systems with switching sequences 
constrained by an automaton. The key element of the framework is the algebraic concept of multinorm, which associates a 
different norm per node of the automaton, and allows to exactly characterize stability. Building upon this tool, we develop the 
first arbitrarily accurate approximation schemes for estimating the constrained joint spectral radius p, that is the exponential 
growth rate of a switching system with constrained switching sequences. More precisely, given a relative accuracy r > 0, the 
algorithms compute an estimate of p within the range [p, (1 -I- r)p\. These algorithms amount to solve a well defined convex 
optimization program with known time-complexity, and whose size depends on the desired relative accuracy r > 0. 
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1 Introduction 

In this paper, we study discrete-time linear switching 
systems having the particularity that their switching se¬ 
quences are constrained by logical rules. We begin with 
an example introducing such systems. 

Given an unstable matrix Ai S and an input-to- 

state matrix B G one computes a control gain 

matrix K G such that A 2 = {Ai + BK) is sta¬ 

ble. The matrix A 2 dictates the closed-loop dynamics of 
a plant, Xt+i = A 2 Xt, whose stability is ensured by a 
state-feedback controller. Let us now consider that the 
controller can fail at any time t, such that the dynamics 
at that time are given by Xt+i = AiXt- Then, the dy¬ 
namics of the plant with failures can be modelled as a 
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switching system 

Xt+l 

where a{t) G {1,2} is the mode of the system and 
(t( 0), cr(l), ..., is the switching sequence that drives the 
system. Without more information on the occurrences 
of the failures, we can only assume that the system is 
unstable. Indeed, in the case of a permanent failure, rep¬ 
resented by the switching sequence a{t) = 1, Vt > 0, the 
plant would follow the unstable dynamics Xt+i = AiXt 
at every time t > 0. However, if we knew with certainty 
that the failure cannot occur more than twice in a row, 
then the above switching sequence would no longer be 
possible, and the system could very well be stable. 

This paper provides tools for the stability analysis of 
switching systems with constrained switching sequences, 
as in the example above. We say that the switching 
system on the matrix set S = {Ai, A 2 ,..., A^} is sta¬ 
ble if and only if, for all accepted switching sequences 
o-Cl)) • ■ we have limt_>oo ^^(t) • ’ ’ ^< 7 ( 0 ) = 0 . 
Switching systems find applications in many theoretical 
and engineering related domains [13,15-17,23,26], and 
the stability of switching systems is known to be a chal¬ 
lenging question [23,24,31]. 

If one does not impose any constraint on switching 
sequences, the resulting system is called an arbitrary 
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switching system. These systems have received a lot of 
attention in the past (e.g., [1,3,4,15]). The stability 
of an arbitrary switching system on a set of matrices 
S is characterized by its joint spectral radius (JSR) 
/5(S) (introduced in [30]). It represents the worst case 
exponential growth rate of the system, and stability is 
equivalent to /5(I]) < 1, which is also equivalent to expo¬ 
nential stability. There has been a lot of research effort 
towards the computation and approximation of the JSR 
(see e.g. [3,8,15,32] and references therein). One com¬ 
mon way to do so is by computing a contractive invari¬ 
ant norm for the system [4,5,8,27], which always exists 
for stable arbitrary switching systems. For any level 
of relative accuracy r > 0, one can approximate these 
norms with quadratic/sum-of-square polynomials [8,27] 
and provide an upper bound on the joint spectral radius 
within the range [p, (1 -I- r)p\. The computation of this 
estimate is done with finite time-complexity. 

Our focus is on the stability of switching systems hav¬ 
ing logical rules on their switching sequences, such as 
the ones studied in [3, 7, 11, 12, 18-21, 24, 29, 33]. We 
refer to these as constrained switching systems, and 
represent the rules by using an automaton. An au¬ 
tomaton is a strongly connected, directed and labelled 
graph G{V,E), with Ny nodes in V and Ne edges 
in E. The edge (v, w,cr) € E between the two nodes 
v,w € V carries the label a G {1, ■ • ■ ,N}, which maps 
to a mode of the switching system. A sequence of modes 
ct( 0), cr(l),..., is accepted by the graph G if there is a 
path in G carrying the sequence as the succession of 
the labels on its edges. We do not specify an initial and 
final node for accepted paths, in that we depart from 
the usual definition for an automaton (see [25], Section 
1.3). The accepted switching sequences form a symbolic 
dynamical system called sofic shift (see [25], Section 
1.5). Examples of automata are given in Figure 1. 



Fig. 1. The labels are represented on the edges. Fig. la corre¬ 
sponds to the example in the introduction, where mode “1” 
cannot occur more than twice in a row. Node “a” is reached 
when the controller works, “b” is reached after one failure, 
and “c” after two failures. The automaton of Fig. lb accepts 
arbitrary switching sequences on N modes. 


The system on the automaton G with matrix set S is 
denoted S'(G, S). The stability of S'(G, S) is character¬ 
ized by the constrained joint spectral radius, introduced 
by Dai [11]. A proof of the following is given in Annex 6. 


Theorem 1.1 (Dai [11], Corollary 2.8) A constrained 
switching system S'(G,S) is stable if and only if its 


constrained joint spectral radius (CJSR), defined as 
p{S) = lim max{llA„(t_i)---A^(o)ll^/* : 

t^oo „{■) (1) 

cr(0),..., a{t — 1) is accepted by G}, 

satisfies p{S) < 1. This also implies exponential stability. 
For all accepted .switching seguences, 


3K > 1,0 < p < 1 : VT > 0, \\A^iT-i) ■ ■ ■ A,,(o) ]] < Kp^. 


The CJSR, defined as (1), is independent of the norm 
used and homogeneous in S. 

To the best of our knowledge, previous works on the 
stability of constrained switching systems have focused 
on establishing algorithmically checkable stability condi¬ 
tions, without studying their conservativeness. There is a 
particular interest in using multiple quadratic Lyapunov 
functions as stability certificates [7,9,12,21,22,24]. These 
approaches provide sets of LMIs whose feasibility is suffi¬ 
cient for stability. In [7,21], a hierarchy of more and more 
complex LMIs is presented such that, for any stable sys¬ 
tem, all LMIs starting from a certain level of complex¬ 
ity (depending on the system) are feasible. The methods 
discussed above can be used to obtain upper bounds on 
the CJSR (a feasible LMI indicates p{S) < 1). However, 
accuracy guarantees on these bounds, similar to that ex¬ 
isting on the JSR estimation, have not been proven yet. 
The framework we introduce allows to obtain accuracy 
guarantees ^ . A direct approach could rely on building 
an arbitrary switching system whose JSR equals the 
CJSR of the constrained system [18,33]. We provide 
more efficient and intuitive techniques. We generalize the 
recent results from [3] towards constrained switching sys¬ 
tems. In [3], the authors focus on systems S'(G, S) with 
G accepting arbitrary switching sequences, and provide 
accuracy bounds for the JSR estimation using multiple 
Lyapunov functions. The generalization of these results 
to general constrained switching systems was left as an 
open question. 

The plan of the paper is as follows. Section 2 introduces 
the algebraic concept of multinorm, which characterizes 
the stability of constrained switching systems as contrac¬ 
tive norms do for arbitrary switching systems. In Sec¬ 
tion 3, we focus on the more algorithmic question of the 
approximation - in finite time and with arbitrary accu¬ 
racy - of the CJSR of a system S'(G, S). In Section 4, 
we illustrate our framework on a numerical example. 
Notations. The matrix A^ G transpose 

of A G A path p of length T > 0 in a graph 

G is a sequence of T consecutive edges. For a path p 
with length T > I, by a slight abuse of notations, we let 
Ap = • • • Acr(i), where the cr(l),..., a{T) are the 

T labels along ». If T = 0, we let A„ = I, the identity 
matrix of R". 


® Preliminary results were presented in [28]. 
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2 Lyapunov functions for constrained switching 
systems 


Proposition 2.2 The constrained joint spectral radius 
(1) of a system S{G, S) satisfies 


The stability of arbitrary switching systems is equiva¬ 
lent to the existence of a contractive norm serving as a 
Lyapunov function. We recall that a norm is a sub ad¬ 
ditive, positive definite and homogeneous function. 

Proposition 2.1 (e.g. [15], Proposition 1.4) The 

joint spectral radius of a set of matrices S is given by 

/5(S) = inf min {y : \Ax\ < 7 |a;|, Va; G R", A G 

|.| 7 t J 

( 2 ) 

where the infimum is taken over all vector norms in R". 


p{S) = inf{ 7 *('H) -.TL is a multinorm for S}. (4) 

PROOF. We first show that the value of any multi¬ 
norm for a system is an upper bound of its CJSR. 
Consider a multinorm TL = G V} for S with 

value 7 . For any path p with length k > 1 between two 
nodes v and w in G, from (3), we get |Apa;|uj < 7 ^|a;|„. 
For any norm j-l, by equivalence of norms in R", there 
exists 0 < a < P such that the inequalities a|a:| < |a;|i, < 
P\x\ hold for all x G R” and all the norms [-It, in TL. 
Considering the classical definition for an induced ma¬ 
trix norm, we have 


A stable arbitrary switching system has /5(S) < 1 (see 
[15]) and from Proposition 2.1, there exists a norm j-j 
such that \Ax\ < [xj, Vx G R", VA G S. It is how¬ 
ever straightforward to build stable constrained switch¬ 
ing systems for which contractive norms do not exist. 

Example 1 A scalar arbitrary .switching system built on 
S = {Ai, A 2 } = {2,1/8} is unstable and it has no con¬ 
tractive norm due to Ai. Consider now the automaton G 
of Figure 2. The periodic system S'(G,S) is stable and, 



Fig. 2. 

applying (1), p{S) = limt^oo ||(AiA 2 )‘||^/(^‘) = 1/2. 

We fill this gap between arbitrary switching and con¬ 
strained switching systems by introducing the algebraic 
concept of multinorm. 

Definition 1 (Multinorm) A multinorm TL for a sys¬ 
tem S'(G(y, E), S) is a set of Ny norms TL = {Hi,, v G 
V}. The value 7 *('H) of a multinorm is defined as 

j^TL) =min{j :\A^x\yj <j\x\y, 

7 (3) 

Vx G R", {v, w, a) G E}. 


Similar ideas have appeared in the literature [2,3,9,10, 
21], where multiple Lyapunov functions are considered 
for characterizing stability where single Lyapunov func¬ 
tions fail to do so. Their role was either to provide a 
sufficient stability condition under the form of a set of 
LMIs, or in [3], to characterize the stability of arbitrary 
switching systems using multiple Lyapunov functions. In 
comparison, we provide general necessary and sufficient 
conditions for the stability of constrained switching sys¬ 
tems, using multiple Lyapunov functions with exactly 
one norm per node of G. 


PpII 


A 


max 

kl=i 



< 


P 

— max 

a \x\=i 


I Apxjuj 
j^ji; 



Taking paths of lengths fc —>■ 00 and the fcth root of the 
above inequality, we obtain p{S) < 7 from (I). 

We now show that for any e > 0 there exists a multinorm 
of value at most {p{S) -\- e). Consider the scaled set of 
matrices 


S' = {A' = A,/(p(5) + e),* = I,...,iV}. 

The CJSR of S'(G,S) is an homogeneous function of 
S (see (I)). The system S'(G,S') is then stable since 
P{S)/{p{S) + e) < 1. 

We define, at each node v G V, the following functions 
which we then prove are norms: 

jxjt, := sup{|ApX| : p is a path with origin v}, 

where j-j is the euclidean norm. These functions are 
well-defined (by exponential stability of ^(G, S')): sub¬ 
additive, homogeneous, and positive definite (with paths 
of length 0, |x|„ > jxj), hence they are norms. Moreover, 
for any edge {v, w, a) G E, and all x G R", we have 

|x|„ = sup{|A/x| : p is a path with origin v}, 

p 

> s\xp{\A'^A'^x\ : <7 is a path with origin ru}, (6) 
9 

= jvlo-xju,, 

where (6) is obtained by taking p starting with the edge 
{v,w,a). Since A],. = A„/{p{S) -\- e), we have [Ao-xjuj < 
{p{S) -\- e)|x|„ for all {v,w,<j) G E. Thus, the value of 
the multinorm with the norms defined in (5) is upper 
bounded by p{S) -\- e. 


We conclude that a constrained switching system is sta¬ 
ble if and only if it has a multiple Lyapunov function 
with exactly one norm per node of its graph. The proof 
of the following is direct from Proposition 2.2. 
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Theorem 2.3 A constrained switching system S{G, S) 
is stable if and only if it admits a multinormH with value 
7 *('H) < 1. Such a multinorm is said to be a Lyapunov 
multinorm for the system. 

Example 2 Consider the switching system S of Exam¬ 
ple 1 with switching constrained by the graph of Figure 2, 
with nodes a and b. Here, the multinorm 

H = {|a;|a, l^lb} = {4\x\, |a:|} , 

where \x\ is the absolute value of x, has a value ofXj^, 
obtained by applying (3) to the system. Thus, the system 
is stable from Theorem 2.3. 

2.1 Extremal multinorms and boundedness 

Given a system S, Proposition 2.2 guarantees that the 
value of any multinorm is an upper bound on the CJSR of 
S. We now investigate the existence of multinorms with 
value egual to the CJSR. The result proposed hereunder 
generalizes the characterization of extremal norms for 
arbitrary switching systems (see e.g [15], Section 2.1.2). 

Theorem 2.4 A system S{G{V, E),!^) admits an ex¬ 
tremal multinorm, i.e. a multinormH* = {[ • |t,,v G V} 
withj*{H*) = p{S), if and only if there exists a constant 
K > 1 such that for any path p in G, 

\Apx\<Kp{Sf-\x\, (7) 

where Tp > 0 is the length of p. 


PROOF. Assume first that p{S) > 0. By homogeneity 
of the CJSR with respect to the set S, we can further 
assume p{S) = 1 without loss of generality (by scaling 
the matrices in S by 1/p{S) > 0). 

We start with the only-if part. Let H* = {HwjU G V} 
be an extremal multinorm. Take any norm j • j. By equiv¬ 
alence of norms in R", there are two scalars a,P > 0 
such that for any node v G V, a|-| < j-jt, < /3|-|. Since 
H is extremal, for any path p between two nodes v and 
w, we have [Apccju, < \x\v from (3), at which point we 
conclude that (7) holds for AT = p/a. 

For the if part, consider a system S with p{S) = 1 and 
define at each node v gV the following norm: 

|a;|* = sup{|Apx| : p is a path in G starting at v}, 

p 

where j • j is e.g. the Euclidean norm. The functions j-j*, 
V G V are indeed norms (|a;|* > |a;| since we also include 
paths of length 0). Also, it is direct to check from the 
equationabove that for any edge (v, ic, cr) G E, lAo-xj)), < 
jxj*. Therefore, the multinorm"H = {j-j*, v G R} is such 
that 7*('H) < 1. From Proposition 2.2, since p{S) = 1, 
we conclude that H is extremal. 


In the case p{S) = 0, we observe that (7) holds if and 
only if Ap = 0 for all path p in G, which is equivalent to 
having the value of any multinorm equal to zero as well. 
Thus, we conclude the proof. 


The result above provides a necessary and sufficient con¬ 
dition for the boundedness of constrained switching sys¬ 
tems with p{S) = 1, i.e. the existence of AT > 1 such that 
for all a; G R" and path p accepted by G, \Apx\ < K\x\. 
The boundedness of systems with p{S) = 1 is known to 
be undecidable (see [29] for sufficient conditions), and 
thus, so is the existence of an extremal multinorm for a 
given switching system. In Subsection 3.4, given a multi¬ 
norm, we provide a sufficient condition for extremality, 
and apply it for computing the CJSR. 

3 Approximation algorithms for stability 
analysis 

Given a system S'(G, S) and a maximum relative error 
r > 0, we wish to obtain an estimate p of the CJSR p{S), 
such that p{S) < p < (1 J- r)p{S). 

Through this section we will provide several approxima¬ 
tion algorithms solving the above problem. All the meth¬ 
ods share the same core mechanism: the approximation 
of a multinorm for the system S, with value close to the 
CJSR, by a quadratic multinorm H = {\ ■ \q,v,v G V}, 
where each norm is quadratic, i.e. [xlg^t, = x^Q^x for 
a positive definite matrix 0. This is expressed as 

a quasi-convex optimization program, solved by using a 
bisection procedure, iteratively checking the feasibility 
of the set of LMIs ( 8 ). 

Theorem 3.1 Consider a system S'(G(R, A), S). 

The value 7 * (A) such that 

7 * (S') = inf 7 

7,{Q„eR"X'“.ijev} 

V {v,w, a) G E, 

-AjQ^A, + ^^Q„>0, ( 8 ) 

Vu G R, ^ 0, 

satisfies the following inequalities: 

KS) < l*iS) < Vn ip{S)). (9) 

Moreover, the LMI feasibility sub-problem (8) is solved 
in a number of operations bounded by 

O , ^ 3 / 2 ^ _ 


PROOF. First, we show that for any system S(G, S) 
and any e > 0 , there is a quadratic multinorm with a 
value 7 that satisfies p{S) < 7 < ^/n{p{S) -\- e). 


4 



The result is obtained from John’s Ellipsoid Theo¬ 
rem [14] (see [8,27] for similar approaches for arbitrary 
switching systems). John’s ellipsoid theorem states 
that for any norm j-j of R”, there exists a quadratic 

norm j-jQ : x —>■ (x^Qx) , with Q 0, such that 
Va; € M” : \x\q < |a;| < ^/n\x\Q. Let us take e > 0 and a 
multinorm He = {Muif G V} with 'j*{He) < p{S) + e. 
Such a multinorm exists (Proposition 2.2). By John’s 
ellipsoid theorem, there exist quadratic norms, forming 
a quadratic multinorm = {Hg,?;,?^ G P}, such that 
for any edge {v, w, a) G S, Vx G M" : 

|^<Ta:|g,u; < \Aax\n, < (/5(S')+e)|a;|„ < v^(p(S')+e)|a;|g,„. 

Since the above holds for any edge we can state that, 
for any e > 0, there is a quadratic multinorm Hq^e with 
7 *(’Hg,e) < yHi{p{S) + e). Taking e —>■ 0 we obtain (9). 
The complexity computations are obtained from the ref¬ 
erence book [6], p.424. The number of variables in the 
problem (9) is (n(n-|- l)/2) The LMIs constraints 
can be represented by a n {Ny -b Ne) bloc diagonal sym¬ 
metric matrix with diagonal blocs of size n x n. 

Notice that 1) the above result gives a sufficient con¬ 
dition for a given system to possess a quadratic Lya¬ 
punov multinorm (if p(S) < Xjsfn^ then (9) guarantees 
7 * (S') < 1) and that 2) it provides an algorithm to solve 
the relative approximation problem with maximum er¬ 
ror r > ^Jn — 1. We now present ways to increase the 
accuracy of the estimation, in Subsection 3.1 through 
Subsection 3.3, by performing an algebraic lifting of the 
structures defining the system S(G, S). 

3.1 The T-product lift 

This method allows to provide arbitrarily accurate esti¬ 
mates of the CJSR, with the cost of constructing a sys¬ 
tem on a graph with a large amount of edges. 

Definition 2 (T-product lift) Given a system S{G,'S) 
and an integer T > 1, the T-Product lift of S, denoted 
S^(G,S), is a constrained switching system on an 
automaton G'{V',E') and a matrix set S' defined as 
follows: 

(1) G' has the same set of nodes as G (i.e. V = V). 
To each path p in G of length T, between two 
nodes v and w in V, is associated an edge 
e = {v,w, {cr(l) ... cr(T)}) G E'. The label on this 
edge is a concatenation of those across the path p. 

(2) The set of matrices S' is the set of all products of size 

T of matrices in S that are accepted byG. For a label 
{cr(l) .. .CT(r)} of the lifted system, = 

Aa{T) ■ ■ ■ 41ct(1) G s'. 

If the system S'(G, S) describes the evolution of a state 
Xt, t = {1, 2,...}, then the system describes the evo¬ 
lution of the same state at times kT for k = {1, 2,...}. 


Example 3 The 2-product lift of the automaton of Fig¬ 
ure 3a is presented on Figure 3b. 

(a-) (b) 

Fig. 3. The graph of Fig. 3b has one edge per path of length 
2 in the graph of Fig. 3a. 

Theorem 3.2 The optimal value^^:{S'^) obtained by ap¬ 
plying Theorem 3.1 to the system S'^{G, S) is such that 
piS) < < n''-A‘^'^)p(S). 


PROOF. We first show that for any system S and any 
integer T > 1, p{S'^) = p{S)'^. Define for all k> 1, 

PkiS) = max{||Tp||^/^ : p of length k accepted by G}. 

p 

Clearly, pk{S'^) = pkriS)'^, and by continuity of the 
exponentiation, we obtain p{S)'^ = limfc_j.oo PkriS)'^ = 
limfc_i.oo /5fc(*S'^) = p(<S'^)- Applying Theorem 3.1 to S'^, 
the inequalities (9) give p{S)’^ = p{S'^) < 7 *(5”^) < 
n^^'^p{S'^) = n^Ap( 5 ')^^ which concludes the proof. 

Corollary 3.3 For any system S'(G, S) and relative 
error bound r > 0, let 

T = riog(n)/(2 log(l -b r))]. 

The CJSR estimation 7 *(S'^) obtained by apply¬ 
ing Theorem 3.1 to the system S'^(G, S) satisfies 
p{S)<j 4S^)^/^ <il + r)p{S). 

As a conclusion, given r > 0, there is a T-product lift of 
S allowing the retrieve an estimate with relative error 
at most 1 -b r. The amount of edges in the system 
increases exponentially with T, with one edge per path 
of length T. The amount of edges Ne of G has however 
the least impact for computational complexity ( 10 ). 

3.2 Approximation through sum-of-squares. 

In this subsection we generalize the technique presented 
in [27] to constrained switching systems and devise a 
CJSR approximation scheme relying on an algebraic lift¬ 
ing of the set of matrices S. The procedure generating 
the lifted space is called the [d]-lift (a full presentation 
can be found in [27], Section 3). Given an integer d, the 
[c?]-lift of a vector x G R" is a vector x'Ai G RC'(n-i-(i-i,d) 
of monomials of degree d, where C{k, £) is the number of 
combinations of i elements in a set of k elements. More 
important to us is that the [d]-lift is well defined for lin¬ 
ear maps, i.e. = (Tx)['^l, this definition 
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extends to sets of matrices Moreover, = |a;|'^ 
for the Euclidean norm in the appropriate dimensions. 
In constrast with the method of Subsection 3.1, we now 
conserve the graph of the constrained system, but ap¬ 
proximate its CJSR by using potentially non-convex ap¬ 
proximations of multinorms that are obtained from ho¬ 
mogeneous sum-of-squares polynomials with degree 2d, 
d > 1 being an integer to be chosen. These polynomi¬ 
als have the form x i—>■ Qx^'^\ where Q 0 is a 

quadratic form in the lifted space. 

Theorem 3.4 Given a system S{G, S) and an integer 
d > 1, let S'[‘^ 1 (G,S) denote the constrained switching 
system on the same graph G and on [(I\-lift o/S. 
The value obtained by applying Theorem 3.1 to 

the system satisfies p{S) < < (j{n + d — 

PROOF. Since holds for the euclidean 

norm, we have = p{SY from (1). Given the di¬ 

mension of the set Theorem 3.1 produces an esti¬ 
mate such that < 7 * < C{n + d — 1, p{S^‘^^). 


The transformation here affects the dimension n of the 
system. As it is shown in [27] - Example 4, the CJSR 
approximation can be further refined by making use of 
explicit sum-of-square programming (rather than solv¬ 
ing the program of Theorem 3.1 in the lifted space), but 
with the same accuracy bounds. 

3.3 Improving accuracy by adding memory to the graph. 


Definition 3 (M-Path-Dependent Lifting) Given 
a system S'(G(V,E),S) and an integer M > 0, the 
M-Path-Dependent lift of S, denoted S'm(G,S), is a 
constrained switching system with the same set of ma¬ 
trices S and with an automaton G'(E', E') which is the 
same as G for M = 0, and constructed as follows for 
M > 1. Start with V = {}, E' = {}, then: 

(1) For each path p of length M in G, add the node Vp 
in V. 

(2) For each path p of length M \ in G, let p = 
(ei,..., Cm+i), where Ck is the kth edge of the path. 
Define p~ = (ei,..., Cm), = (e 2 ,..., Cm+i), 
and Op as the label of the edge Cm+i- Then, add the 
edge ivp-,Vp+,ap) G E'. 

Example 4 Figure 4 presents the 1-Path-Dependent lift 
G' of the automaton G of Figure 3a. There are 4 nodes 
in G', one for each edge ofG. For the edge {v,w,a) of 
G, we use “vw” to refer to the corresponding node in G'. 
For the edges, consider for example the path of length 2 
((6, b, 2), (5, a, 4)) inG. To this path corresponds the edge 
{bb, ba, 4) in G'. 



Theorem 3.5 Given a system S'(G,S) and an in¬ 
teger T > 1, let 7 *(S'^) and 7 *(S't-i) be the val¬ 
ues obtained by applying Theorem 3.1 to S'^(G, S) 
and S't-i(G,S) respectively. The estimates satisfy 


Path-dependent Lyapunov functions have been intro¬ 
duced by Lee and Dullerud [21] as tools for the stability 
analysis and control of discrete-time switching systems. 
The concept, which follows a similar idea to that of Bli- 
man and Ferrari-Trecate [7], is to build a multiple Lya¬ 
punov function that associates a different quadratic form 
to each switching sequence of a length M > 0, M being 
an integer parameter called the memory of the function. 
The authors showed that, for any stable switching sys¬ 
tem, there is a finite M such that the system admits a 
path-dependent Lyapunov function with memory M. 
Similar statements can be made about the T-Product 
lift and the [dj-lift defined in the previous subsections. 
Indeed, if p{S) < 1 and since the approximations are 
asymptotically tight, there is a finite value T (or d) for 
which the approximation algorithm will return an esti¬ 
mate lower then 1. Given this parallel, it is natural to 
ask whether the methods of [7, 21] present similar ap¬ 
proximation bounds to that of Theorems 3.2 and 3.4. 
We first define a lifting procedure, adapted from [7,21], 
allowing to obtain path-dependent Lyapunov function 
as quadratic multinorms of an augmented automaton. 


PROOF. For T = 1, S' = St-i = S^, so the claim 
holds. Assume now T > 2. We will show that given a 
quadratic multinorm for S^ with value 7 ^, we can al¬ 
ways construct a quadratic multinorm for St-i with 
value at most 7 . 

We refer to paths by their succession of edges, p = 
(ei,..., ct) is a path of length T in G, and let Ck = 
{vk-i,Vk, CTfe). A quadratic multinorm with value 7 ^ for 
S^ associates a quadratic form Q >~ 0 per node in G 
such that, for all paths p, 

( 11 ) 


To the same path p corresponds an edge in the (T-1)- 
Path-Dependent lift, between a node V(^ei,...,eT-i) ^ 

node F'or a quadratic multinorm of St-i with 

value 7 to exist, there must be a quadratic form Qp 0 
per path p of length T, such that the following holds: 


k^(TT^(e2,...,eT)^f7 




( 12 ) 
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Letting R = Q~^ denote the inverses of the quadratic 
forms in the corresponding lift, by a Schurr complement, 
the LMIs (11) are equivalent to 


^pRvQ^p q Rvt 1 — O5 

and the LMIs (12) to 

^( 7 t “^(ei ,...,eT —1 )^<7t 

^ 0 . 


(13) 

(14) 


Assume that we have a solution {Rv,v G I^} to 
(13). Given any path p = (ei, ..., Ct-i), with labels 
CTi,..., (Tt-i that visits the nodes vq, ■ ■ ■, vt-i in G, 
define 


D _ D 

~ ^ vt -1 

+ 7 ‘^^aT-l^VT-2^T-l 

7 ^((Tt-2,'7T-i)'^1'T-3 


A 


T 


+ 7 


where = ^^t-i ’ ’ ’ ■ Injecting these 

quadratic forms in (14), we obtain (13), and thus we 
have a quadratic multinorm for St-i, with value at 
most 7 . 


Corollary 3.6 Given a system S'(G, S) and an integer 
M > 1, the value 7 * obtained by applying Theorem 3.1 
to the M-Path-Dependent lift of S satisfies p{S) < 7 * < 
ni/( 2 (^+i))p(^). 


The amount of nodes Ny and edges Ne in the M-Path- 
Dependent lift both grow exponentially with M, with 
one node per path of length M, and one edge per path 
of length M + 1 in G. 


3.4 A sufficient condition for extremality of quadratic 
multinorms. 


We now present an easy to check sufficient condition un¬ 
der which the CJSR estimate obtained by applying The¬ 
orem 3.1 is exact. The condition can be applied for all 
the lifts developed above, since they all rely on Theorem 
3.1 to retrieve a quadratic multinorm of minimal value. 
We assume that the estimate of Theorem 3.1 is attained 
by the value of a quadratic multinorm, and give a suf¬ 
ficient condition for its extremality. We start with the 
following observation. 

Lemma 3.7 For any system ^(G, S) and any cycle'^ c 
of length T in G, p{S) > p{AcY/'^, where p{Ac) is the 
spectral radius of the product Ac. 

^ A cycle is a closed path, i.e. whose source and destination 
nodes are the same. 


PROOF. For any induced matrix norm || • ||, the fol¬ 
lowing holds 

p(A,)i/^= hm WAlW^/^^ 

k—¥OQ 

< lim max{|| : p is a path of length fcT}, 

k—¥oo 

< KS)- 


We now define a simple cycle in the graph G as a cycle 
such that for every node v visited by the cycle, there 
exists no partition of the cycle into two cycles on v. 

Theorem 3.8 (Sufficient extremality condition) 

Let{Qv, V G V} be the optimal quadratic forms obtained 
by applying Theorem 3.1 to the system S{G(V, E),!!), 
corresponding to a multinorm TL with value 7 *. If the set 
of edges 

E' = {{v,w,a) & E : A^in {-AJQ^A^-\-j‘^Q^) = O} , 

where Amin(-A) denotes the smallest eigenvalue of the 
positive definite matrix X, forms a simple cycle c in G, 
then TL is extremal, and p{S) = p(Ac)^/^, where T is 
the length of c. 


PROOF. The edges in the set E' are those for which 
the LMI constraints in ( 8 ) are tight for the multinorm 
"H. It is a known fact of convex optimization that remov¬ 
ing constraints that are not tight at the optimal solution 
of a given program does not affect its optimal solutions. 
Therefore, it must be the case that the optimal values 7 * 
and 7 ( obtained by applying Theorem 3.1 respectively 
to S'(G(F, E),!^) and S{G{V, A'), S) are equal, and we 
may focus on the second system. 

Since G{V,E') is a simple cycle, then for any M > 0, the 
M-Path-Dependent lift (Definition 3) leaves the graph 
invariant. Applying Corollary 3.6 with M —>■ 00 , we con¬ 
clude that the multinorm obtained from applying Theo¬ 
rem 3.1 to S'(G(F, A'), S) needs to be extremal for this 
cyclic graph, and 7 ^ = p(Ac)^/^. 

Having computed 7 ^,, and since 7 * = 7 ^, we now consider 
the original system. Applying Lemma 4 and Proposition 
2.2, we get p{Ac)^^'^ < p{G, S) < 7 * = p(Ac)^/^, which 
concludes the proof. 


3.5 The Kroenecker lift 

We end this section by presenting another approach to 
the stability analysis of constrained switching systems. 
Kozyakin [18] and Wang [33] independently introduced a 
lifting procedure that allows to obtain, from any system 
S, a set of matrices Eg such that p{S) = p{T^s)- 


7 



Definition 4 (Kroenecker lift) Given a system 
S{G{V, with Ny nodes {vi & V, 1 < i < Ny} 

and S C the Kroenecker lift of S is a matrix set 

defined as Sg = {vi,Vj,a) € E}, with for 

the edge {vi,Vj,a), 

where e{k) G is the kth vector of the canonieal basis 
o/R^'^, and 0 is the Kroenecker product. 

The methods in Subsection 3.1, 3.2 and 3.3 are based 
on the ones of [8,21,27] and make use of the concept of 
multinorm for providing estimates of p{S). It is natural 
to compare these methods with a direct application of 
the ones of [8,21,27] for approximating /5 (Sg). We will 
focus on [ 8 ], which approximates the JSR of a set of 
matrices using Theorem 3.1 with a single quadratic form. 
The conclusions naturally carry on to the other methods 
of [21,27]. 


using Sg C ^ complexity ( 10 ) reads 

a significant increase compared 

to the previous O . Moreover, the 

accuracy bounds (9) become worse, due to the larger 
dimension of Sg- 

4 Numerical example 

We consider the dynamics of a plant that may experience 
controller failures: Xt+i = (^ + xt, with 

0.94 0.56\ 

0.14 0.46/ Vl/’ 

and K„(^t) = The control gains switch 

to represent 4 failures modes. When everything works as 
expected, aft) = 1 , and 


Proposition 3.9 Consider a system S{GfV,E),T,) 
and its associated set Sg- There exists Qg >- 0 and 
7 > 0 such that'i G Sg, 

T Qg 0, 


if and only if S has a quadratic multinorm with value at 
most 7 . 

PROOF. Assume without loss of generality that 7 = 
1. To ease the reading, we assume n = 1, so that for 
e G E, Af. G R-^vxAfy^ Pqj. g _ yfQ let 

Be = AjQf^Ae — Qg- The element at row k and column 

i of Ae is and bI!^’^\ are defined similarly. 

Define the set of nodes as G F, 1 < f < IVy}. For 
any I < k,£ < Ny, we have 


(vi,Vj,cr) 




T,[fc,i'] fAi',3'\ aM'A 


i'd' 


^[k,£] 



T.[fcj] ^[jj] 

{vi,vya)AG 


{Vi,Vj,(7) 


Vg ’ 


where Sij = 1 if and only if i = j (Kroenecker delta). 
For the only if part, Qg ^ 0 and Be E 0. Since Qg ^ 0 
and Be E 0, then for al\l < k < Ny^ Qq > 0 , and for 

any e = y,Vj,a) G E, M*’*' = AjQ^y^Aa^ - < 0. 

Thus, we extract a quadratic multinorm for S using the 
quadratic forms Qg*^, 1 < i < Ny 

For the if part of the proof, we take quadratic multinorm 
with value lower than 1. Let Qi > 0,1 < f < Ny, be the 
quadratic forms associated at each node. We reconstruct 
Qg by setting Qq = Q;. if fc = ^, and Qq = 0 else. 


Ki = (/ci fe) = (-0.49 0 . 27 ) . 

The second and third modes correspond respectively to 
a failure of the first and second part of the controller, 
with K 2 = ( 0 ,^ 2 ) and K 3 = (A:i,0). The last mode is 
the total failure case with K 4 = (0,0). As a constraint, 
we consider that a same part of the controller never fails 
twice in a row. We obtain a constrained switching sys¬ 
tem on the set S = {A -|- BK^}, a G {1,..., 4}, with 
the automaton G(F, E) depicted in Figure 5. 

We first provide increasingly accurate estimates of the 



Fig. 5. 

CJSR of the system, and then provide the exact value 
of the CJSR using Theorem 3.8. For the estimations, we 
compare the T-Product lift of Subsection 3.1, and the 
M-Path-Dependent lift, introduced in [21] and further 
studied in Subsection 3.3. We apply these methods for 
T = 1,..., 7, and M = T — 1 = 0,..., 6 . From Corol¬ 
lary 3.3 and Theorem 3.6, these choices produce esti¬ 
mates within [/5(S'), (1 -|-■^’)/5(S')], with 1 -|-r ranging from 
V2 ~ 1.41 for T = 1 to 21/14 _ Q 5 p = 7 , 

For each value of T and M we first compute the lifted 
systems and Sm- Figure 6 provides a comparison of 
the amount of nodes and edges of these systems. Then, 
we solve the optimization program of Theorem 3.1 on the 
lifts, applying Theorems 3.2 and 3.5 to obtain bounds 
on the CJSR. ® The execution times for producing the 


Applying Theorem 3.1 to Eg therefore produces the ^ Matlab codes reproducing the results available at 
same estimate as if it was applied on S. However, when “http://sites.uclouvain.be/scsse/postrev.zip”. 



estimates are compared in Figure 7, and the estimation 
values are compared in Figure 8. The best estimation is 
achieved with the M-Path-Dependent lift, for M = 6, 

guaranteeing that p{S) < 0.9748-However, the T- 

Product lift provides estimates much faster. 

We now use Theorem 3.8 for computing the exact 




Fig. 7. 



Fig. 8. Using the M-Path-Dependent lifts, we produce better 
estimates. The estimates produced by the T-Product lifts 
need not decrease monotonically as T increase. Lower bounds 
computed from the relative accuracies of the CJSR estimates. 

value of the CJSR. The results of Figure 8 indicate that 
the sufficient condition can only be met by the M-Path- 
Dependent lift, for M = T — 1 = 4,5 or 6. In this 
case, for M = 5, the conditions of Theorem 3.8 are met. 
The obtained simple cycle corresponds to the switch¬ 
ing sequence repeating the labels {2, 3,1,1,1,1, 2,1} and 


reaches the CJSR. Indeed, we have (Lemma 4) 

KS) > piAiA2AjA3A2y/^ ~ 0.9478..., 
which matches our best CJSR estimate. 


5 Conclusion 

We established a new framework for the stability 
analysis of discrete-time linear switching systems with 
constrained switching sequences. It relies around the 
newly introduced concept of multinorms and their link 
to the constrained joint spectral radius [II] (CJSR). 
The stability of constrained switching systems is equiv¬ 
alent to the existence of a set of exactly one norm per 
node of the automaton, with contractivity relations that 
are given by the edges of the automaton. 

By approximating these norms individually with 
quadratic norms, we design the first arbitrarily accurate 
CJSR approximation schemes. The framework also en¬ 
capsulate well-known methods such as path-dependent 
Lyapunov functions [7, 21] as methods approximating 
extremal multinorms. 

In the future, we will apply this framework to control- 
oriented problems. Path-dependent Lyapunov functions 
have been used in feedback controller design [12,21], and 
our goal is now to give guarantees on the performance 
of these controllers. Their usage for the study of the 
contractiveness of switching systems [12,20] suggests an 
estimation framework in this setting. 
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6 Annex: Proof of Theorem 1.1 

From Definition 2 in [11], we know that the limit (1) 
converges. Since the limit (1) converges, Ve > 0, > 

0 : Vt > T„ 

KSy < max ||A<,(i_i)---A,^(o)|| < (/5(S')+e)*. 

(7(-) accepted by O 

To show that p{S) < 1 implies (exponential) stability, it 
suffice to take e < 1 — p{S). Indeed, we then get that for 
all accepted sequences, limt^oo ll^o-(i-i) • • • ^<t(o) II = 0- 
Exponential stability is then acquired since there can 
only be a finite amount of products of length t < Te with 

||clcr(t-l) • • • II > (p('S') + e)*- 

Consider now the case p{S) > 1. For alH > 1, we define 
x\ as 

x*=argmax max \A„(t-i) ■ ■ ■ 

|a:| = l (7(-) accepted by G 

We extract from the sequence {a;‘}t>i a subse¬ 
quence converging to a point x*, |x*| = 1. From 
this point, there is a switching sequence satisfying 
limt^.oo l^o-(t-i) • ■ ■ Acr{o)X^,\ > 1/2. Thus, the system is 
not asymptotically stable when p{S) > 1. 

This concludes the proof of Theorem 1.1. 
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